Violation of Wiedemann-Franz law at the Kondo breakdown quantum critical point 
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We study both the electrical and thermal transport near the heavy-fermion quantum critical point 
(QCP), identified with the breakdown of the Kondo effect as an orbital selective Mott transition. 
We show that the contribution to the electrical conductivity comes mainly from conduction elec- 
trons while the thermal conductivity is given by both conduction electrons and localized fermions 
(spinous), scattered with dynamical exponent z = ?>. This scattering mechanism gives rise to a 
quasi-linear temperature dependence of the electrical and thermal resistivity. The characteristic 
feature of the Kondo breakdown scenario turns out to be emergence of additional entropy carriers, 
that is, spinon excitations. As a result, we find that the Wiedemann-Franz ratio should be larger 
than the standard value, a fact which enables to differentiate the Kondo breakdown scenario from 
the Hertz-Moriya-Millis framework. 
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Existence of quasiparticles is the cornerstone in Lan- 
dau's Fermi-liquid theory [l| for modern theory of met- 
als. Since they transport not only electric charge but 
also entropy, one sees that the ratio (L = ^) between 
thermal (k) and electrical (ct) conductivities is given by a 

universal number Lq ^ ^(^)^ 2.45 x V^-^W^R-"^ 

provided quasiparticles do not lose their energy dur- 
ing collisions, and certainly satisfied at zero temperature 
in the Landau Fermi liquid theory. Not only conven- 
tional metals [s'l but also strongly correlated metals such 
as heavy fermionsM turn out to follow the Wiedemann- 
Franz (WF) law [5|. In particular, even quantum critical 
metals oi CeN^Ge-i, 0, CeRhln^ 0, and CeCoIn^ 
have shown to satisfy the WF law at least in the zero tem- 
perature limit, thus validating the quasiparticle picture, 
although their resistivities deviate from the conventional 

behavior. 

Several years ago, violation of the WF law was 
observed in the optimally electron-doped cuprate 
(Pr, Ce)2CuO/^ Q] and hole-underdoped cuprate 
La2-xSrxCu04 |9I| while the WF law turns out to hold 
in the overdoped cuprate Tl2Ba2CuOQ+s suggest- 
ing emergence of non-Fermi liquid physics as proximity 
of a Mott insulator. Recently, anisotropic violation of 
the WF law has been reported near the quantum critical 
point (QCP) of a typical heavy-fermion compound 
CeCoIn^, where only c-axis transport violates the WF 
law while ab-plane transport follows it pT| . In this 
experiment the authors speculated that the temperature 
quasi-linear electrical resistivity and vanishing spectral 
weight may be one common feature for such non-Fermi 
liquid physics. In this respect they naturally proposed to 
observe violation of the WF law in YbRh2Si2 as another 
typical hcavy-fcrmion compound exhibiting similar 
non-Fermi liquid physics with the c-axis measurement 
of CeCoIn^, where both ab- and c- axis transport show 
the temperature quasi-linear electrical resistivity and 
vanishing spectral weight p^ . 



Physically, one can expect violation of the WF law 
as proximity of Mott physics or superconductivity away 
from quantum criticality, and as emergence of non-Fermi 
liquid physics near QCPs. In a Mott insulator the pres- 
ence of charge gap makes electrical conductivity vanish, 
but gapless spin excitations can carry entropy, causing 
L > Lq while Cooper pairs transport electric currents 
without entropy in the superconducting state, resulting 
in L < Lq. On the other hand, entropy is enhanced near 
QCPs due to critical fluctuations, and violation of the 
WF law is expected in principle. 

In this paper we examine thermal transport and vi- 
olation of the WF law based on the Kondo breakdown 
scenario [l^, as one possible heavy-fermion quantum 
transition for YbRh2Si2- This scenario differs from the 
standard model of quantum criticality in a metallic sys - 
tem, referred as the Hertz-Moriya-Millis framework p^ . 
in respect that in the former case the whole heavy Fermi 
surface is destabilized at the QCP. 

Several heavy-fermion compounds have been shown 
not to follow the SDW theory [H, [H, [11. Strong 
divergence of the effective mass near the QCP [ll and 
the presence of localized mag netic moments at the transi- 
tion towards magnetism [l7| seem to support a more ex- 
otic scenario. In addition, rather large entropy and small 
magnetic moments in the antiferromagnetic phase may 
be associated with antiferromagnetism out of a spin liq- 
uid Mott insulator [l^ . Combined with the Fermi surface 
reconstruction at the QCP [H, this quantum transi- 
tion is assumed to show breakdown of the Kondo effect as 
an orbital selective Mott transition [l^ , where only 
the f-electrons experience the metal-insulator transition. 

Our main result is that the WF law should be violated 
at the Kondo breakdown QCP as proximity of spin liquid 
Mott physics, thus L > Lo, resulting from the presence of 
additional entropy carriers, here spinon excitations. This 
result may be opposed to the preservation of the WF law 
in the SDW framework 
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We start from the U(l) slave-boson representation of 
the Anderson lattice model (ALM) in the large-C/ limit 

Lalm = c]^{dr - ^J)cia - i ^(cLcj<T + H.c.) 

i (ij) 

+ VJ2ib^flc^a + H.C.) + ^ bldrb, 

i i 

+ E fl(^r + 'f)f^- + jT.^flXvfj<r + H.C.) 

+^ E ^Mb. + flha - 1) + iVJ ^ Ix., r (1) 

Here, Ci^ and di^ = ^I/io- ^^'^ conduction electron with a 
chemical potential and localized electron with an en- 
ergy level e/, where bi and fia- are holon and spinon, 
associated with hybridization and spin fluctuations, re- 
spectively. The spin-exchange term for the localized or- 
bital is introduced for competition with the hybridization 
term, and decomposed via exchange hopping processes of 
spinons, where Xij is a hopping parameter for the decom- 
position. Xi is a Lagrange multiplier field to impose the 
single occupancy constraint bjbi + fj^fia = N/2, where 
N is the number of fermion flavors with cr = 1, A^. 

Performing the saddle-point approximation of bi — > b, 
Xij ~^ Xj and iXi — > A, one finds an orbital selective Mott 
transition as Kondo breakdown at J w Tk, where a spin- 
liquid Mott insulator {{bi) = 0) arises in J > Tk while a 
heavy- fermion Fermi liquid {{bi) ^ 0) results in Tk > J 

Here, Tk = D e^p(^j^jfy^^ is the single-ion 

Kondo temperature, where pc ~ {2D)^^ is the density of 
states for conduction electrons with the half bandwidth 
D. 

One can read the WF ratio in the mean-field approxi- 
mation. In the heavy-fermion phase it is given by L = Lq, 
representing a Fermi liquid state of heavy quasiparticles 
with a large Fermi surface. On the other hand, it be- 
comes 



t + JxV 



t 



in the spin liquid phase, where t is the hopping of the 
conduction electrons and x is the spin liquid parameter. 
By contrast, in the U(l) slave-boson mean- field theory of 
the t-J Hamiltonian, 



Hir 



Y,{NJ\X^,\' - its + JXr,)flU - H.C.y 



where the holon field is replaced with its mean-field 
value of bi = ^/S with hole concentration 6, one finds 

L = Lq ( ) \2l\ . which represents a strong viola- 



te \^ Na- 
tion of the WF law at the vicinity of the insulating phase. 
This comparison tells us that the orbital selective Mott 
transition in the ALM has milder violation of the WF 
law than the single-band Mott transition although the 



underdoped state of the t-J model may have similarity 
with the fractionalized Fermi liquid UJj of the ALM. 

Fluctuation-corrections are treated in the Eliashberg 
framework T^. The main physics is that the Kondo 
breakdown QCP is multi-scale. The dynamics of the hy- 
bridization fluctuations is described hy z = 3 critical the- 
ory due to Landau damping of electron-spinon polariza- 
tion above an intrinsic energy scale E* , while by z = 2 
dilute Bose gas model below E*, where z is the dynam- 
ical exponent. The energy scale E* originates from the 
mismatch of the Fermi surfaces of the conduction elec- 
trons and spinons, shown to vary from 0(10°) mK to 
O(IO^) mK. Based on the z — 3 quantum criticality, a 
recent study [ll] has fitted the divergent Griineisen ratio 
with an anomalous exponent 0.7. 

Transport coefficients can be found from the following 
transport equations 



T 



(2) 



"^eilth) electric (thermal) current for conduction elec- 
trons, spinons, and holons, respectively, and E, e, Hcj.b, 
and T are an external electric field, internal one, each 
chemical potential, and temperature, respectively, where 

and 



ac,f,b = 1,0,-1 and f3c.j,b = 0,1,1. K, 



K2^"^ are associated with electrical conductivity, thermo- 
electric conductivity, and thermal conductivity for each 
excitation, respectively. Obtaining e from the current 



constraint J^^ -I- J^; = with /^c = fj,f — fib, and consider- 
ing the open-circuit boundary condition, we find physical 
response functions for electrical conductivity at, thermo- 
electric conductivity pt, and thermal conductivity Kt, 



jb 



Ct = CTc 



^b(Jf 
(7b + CTf' 



Pt^Pc 



Kt_Kc K/ Kfc_ {Pb+Pf)^ 

r~r T r ab + aj 



(JbPf - CTfPb 
CTb+CTf 



(3) 



with (7, 



cj,b 



PcJ,b 



K{^' /T, and Kc,f,i 



K2 /T. One can also derive this general expression 
from the path-integral representation with the covariant 
derivative D = V — lA^i — iAth{idr) in the continuum ap- 
proximation, where A^i and Ath are external electromag- 
netic and thermal vector potential fields, respectively. 
We note that this expression reduces to that of the t-J 
model, when contributions from the conduction electrons 
are neglected [23j . 

It is straightforward to evaluate all current-current cor- 
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relation functions in the one loop approximation. We find 



a,{T)^Cp,vfTlJT), 



TT 



T 



TT 



(4) 



with C = ^ I^oo '^y (y'^+i)'^ • electrical conduc- 

tivity Pc{f) and Vp'^^ are density of states and Fermi 
velocity for conduction electrons (spinous), respectively. 
'''c(f) scC^) ~ [^'^c(f){T)]~^ is the scattering time due to 

yb 

Qs,(^)(r>£;*) 



z = 3 hybridization fluctuations, given by 

r2 



/(c)- 



E*J 



12t:v 



fie) E* 



In 2, 



{2NV'^Pc) ^ is the band mass for holons. 



3 1 is the transport time as- 



where mi, 

Note that hybridization fluctuations are gapped at 
T < E* , resulting in the Fermi liquid like correction. 

sociated with z = 3 gauge fluctuations, where -fa ~ '^/''^p 
is the Landau damping coefficient for gauge fluctuations 
and kp is the Fermi momentum of spinous. In the ther- 
moelectric coefficient ep is the Fermi energy for con- 
duction electrons, and cp is a geometrical factor, here 
cp = 3/2 for the spherical Fermi surface [24|. 

Several remarks are in order for each transport coeffi- 
cient. An important point is that the vertex corrections 
for scattering with hybridization fluctuations can be ne- 
glected, a unique feature of the two band model, resulting 
from heavy mass of spinous [H, ■ This allows us to 
replace the transport time with the scattering time for 
such a process. On the other hand, vertex corrections for 
scattering with gauge fluctuations turn out to be crucial, 
where infrared divergence of the self-energy correction at 
finite temperatures is cancelled via the vertex correction, 
giving rise to gauge-invariant [1^ finite physical conduc- 
tivity [2^ . As a result, the gauge non-invariant diver- 
gent spinon self-energy 9I]^(T) in 9i;^(T) + '^sY.'jiT) of 
the conductivity expression is replaced with the gauge 
invariant finite transport time [tj f.^(T)]~^ . Both irrele- 
vance (hybridization fluctuations) and relevance (gauge 
fluctuations) of vertex corrections can be also checked in 
the quantum Boltzman equation study. 

Both the thermoelectric and thermal conductivities are 
nothing but the Fermi liquid expressions, where each 
fermion sector satisfies the WF law. Although inelas- 
tic scattering with both hybridization and gauge fluctu- 
ations may modify the Fermi liquid expressions beyond 
the one-loop approximation, the WF law for each fermion 
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FIG. 1: (Color online) Left: Electrical conductivity from con- 
duction electrons (blue), spinous (red), and holons (green). 
Left-inset: Electrical conductivity from holons much smaller 
than contributions from fermions. Right: Thermoelectric 
conductivity from conduction electrons (blue), spinous (red), 
and holons (green). Right-inset: Thermoelectric conductivity 
from holons much smaller than contributions from fermions. 



sector will be preserved at least in the zero temperature 
limit, where such inelastic scattering processes are sup- 
pressed. One may regard the WF law for each fermion 
sector as the most important assumption in this paper. 

Transport coefficients for holon excitations turn out 
to be much smaller than fermion contributions, that is, 
iciT) > af{T) » ab{T), p,{T) > pf{T) » pfc(T), and 
i^c{T) > Kf{T) 3> K}j{T) as clearly shown in Fig. 1, thus 
irrelevant. Physically the dominance of fermion contri- 
butions can be understood from an argument of density 
of states. Since there are many states at the Fermi sur- 
face in the vacuum state, their conductivities diverge in 
the clean limit as the temperature goes down to zero. On 
the other hand, there are no bosons at zero temperature, 
thus their conductivity vanishes when T — > 0. 

Inserting all contributions into Eq. (3), we find the 
physical transport coefficients near the Kondo breakdown 
QCP. Interestingly, the dominance of fermion contribu- 
tions allows us to simplify the total transport coefficients 
as 



at{T)^a,{T)=Cp,vfT-,AT), 
6 ep 



T 



T 



+ 



T 



^(a,(r)+a^(T)).(5) 



Actually, we have checked that each approximate formula 
matches with each total expression. The main point is 
that spinous participate in carrying entropy, enhancing 
the thermal conductivity, while both electric and thermo- 
electric conductivities result from conduction electrons 
dominantly. 

Fig. 2 shows the quasi-linear behavior in tempera- 
ture for both electrical and thermal resistivities above 
E* , resulting from scattering with z = 3 hybridiza- 
tion fluctuations dominantly, because of [t^ g^{T)]^^ ^ 
.(T)]~^ in the spinon conductivity, explicitly checked 



from numerical analysis and temperature dependence, 
thus (Jf{T) w Cpfv{^Tf ,^(T). Here, we have used pa- 
rameters of Ref. [22], shown to be successful for fltting 
of Griineisen ratio. 
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FIG. 2: (Color online) Quasi-linear electrical (blue) and ther- 
mal (red) resistivity above E* , where thermal resistivity is 
smaller than electrical resistivity owing to the contribution 
from spinon excitations. 



The WF ratio is given by 

Kt{T) Kc{T)- 



L{T) 



Tat{T) 



Tac{T) 



PcVp 



(6) 



in the low temperature limit, where gauge- fluctuation 
corrections are irrelevant compared with hybridization- 
fluctuation corrections, thus ignored in the last expres- 
sion. 

The larger value of the WF ratio is the characteristic 
feature of the Kondo breakdown scenario, resulting from 
additional entropy carriers, here the spinon excitations. 
If we perform the transport study based on the SDW 
theory in the same approximation as the present frame- 
work, we win find atiT) cc Ttr{T), ptiT) = ^^Tut{T), 



and 



«t(T) 



^fTt(T), where likewise contributions from 



critical boson excitations are assumed to be irrelevant, 



and the scattering time is replaced with the transport 
time. As a result, the WF law is expected to hold al- 
though non-Fermi liquid physics governs the quantum 
critical regime. Actually, this has been clearly demon- 
strated in the self-consistent rcnormalization theory, well 
applicable to CeNi2Ge2 6]. In this respect the viola- 
tion of the WF law discriminates the Kondo breakdown 
scenario from the SDW framework. 

In this study we have ignored thermal currents driven 
by gauge fluctuations, known to be the phonon-drag ef- 
fect in the electron-phonon system [23|. Recently, Nave 
and Lee have considered such photon-drag effects in the 
spin liquid context with z = 3 gauge fluctuations based 
on quantum Boltzmann equations, and argued that drag 
effects are subdominant, compared with fermion contri- 
butions |27j| . Resorting to their conclusion, we argue that 
dominant thermal currents are driven by fermion excita- 
tions, here both conduction electrons and spinous. 

In conclusion, we found marginal Fermi liquid physics 
for both electrical and thermal transport near the Kondo 
breakdown QCP due to scattering with z = 3 hybridiza- 
tion fluctuations. Our main discovery is that the Kondo 
breakdown QCP should violate the WF law at least in 
the zero temperature limit due to proximity of spin liquid 
Mott physics, i.e., existence of additional entropy carri- 
ers, that is, spinous. 
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